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Some recent publications by authors from the University of Maryland analyse the fluctuations between multi- 
port model parameters in stochastic environments. The authors use random matrix theory (RMT) for estimates 
concerning eigenfunction expansions in cavities which show so-called "wave chaos." A specific and "universal" 
relation between the variances of the fluctuations of multi-port impedance parameters, of the same type as we 
had observed experimentally for scattering parameters in mode stirred chambers, is analysed in detail. However, 
contrary to our own observations these recent papers claim that only the impedance fluctuations satisfy universal 
variance ratios. In this paper, we detail our own approach to the problem which is based on electromagnetic 
scattering theory. It is shown that universal variance ratios are obtained, for any type of multi-port model, when 
the environment of the multi-port system is described by a scattering operator, which is statistically isotropic in 
a sense we make precise below. The only condition on the system is that the radiation power correlation between 
the ports considered vanishes. This is to be compared to the diagonality of the radiation impedance required for 
the Maryland results to be valid. 
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I. INTRODUCTION 

The theoretical research we present in this paper has been 
inspired by some remarkable experimental evidence showing 
the existence of a very specific correlation between the fluctu- 
ations of multi-port S-parameters in a Mode Stirred Chamber 
(MSC). When we were trying to establish precise statistics of 
antenna response in MSC, it appeared that the standard mea- 
surement equipment of an MSC did not allow us to achieve 
the desired accuracy. So the decision was made to change 
the standard measurement setup, with generators, power am- 
plifiers and spectrum analysers, for a network analyser based 
system. In this way, we could indeed establish detailed and 
accurate statistics, over a large number of mode stirrer ori- 
entations, of the complex two-port S-parameters, describing 
the interaction between an emitting antenna and a receiving 
antenna. When plotting the variance curves as a function of 
frequency, it appeared that the variances of the transmission 
coefficients were related to those of the reflection coefficients 
in a very special way. 

The mathematical expression of the experimentally ob- 
served relation was rapidly found to be 

var(5i 2 ) = ^v/var(Sn)var(S 2 2) (1) 

Figure Q] shows the experimental curves as well as the curve 
of var(5i2) computed from the var(Sn) and var(S22) using 
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the above equation. It can be observed that the computed 
curve is almost indistinguishable from the measured one for 
frequencies above 6 GHz. These relations have afterwards 
been confirmed by similar measurements in other MSC and 
mode-stirred enclosures and for various different emitting and 
receiving systems. 

These observations have been published in 10] and the es- 
sentials of the electromagnetic foundation of these relations 
have beenpresented in the Ph.D. thesis of C. Fiachetti [2] 
(see also J5[], which contains another theoretical explanation). 
However, some recent work by researchers from the Univer- 
sity of Maryland ll 111 , has revived our interest in the subject, 
in particular because these authors apparently come to conclu- 
sions contradicting our results. 

If we try to understand the relation between S-parameter 
statistics, the first thing to observe is, obviously, that the var- 
ious 5-parameters appearing in the coupled antenna configu- 
ration must be correlated in some way or another. After all, 
energy which is not emitted, due to an impedance mismatch 
on the emitting side, cannot arrive at the receiving side. Simi- 
larly, the fraction of the incident energy capted by the receiv- 
ing antenna will be affected by the mismatch on the receiving 
side. First attempts to derive the relation in equation ([T]i as a 
simple consequence of energy conservation and/or reciprocity 
failed, though. As it appears, the field-theoretical derivation 
requires a more profound analysis of the electromagnetic in- 
teraction between linear multi-port systems and a stochastic 
environment (see 10, 0] for more background information). It 
is the purpose of this paper to show how this analysis can be 
done. 
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FIG. 1: Experimental evidence for a hard relation between the vari- 
ances of multi-port S-parameters in MSC. The dashed curve is com- 
puted from the two reflection parameters using eq. (TJ. 



Outline of the paper 

In section|nl we summarise the essentials of frequency do- 
main electromagnetic theory as far as we need it in our analy- 
sis. In particular, we formulate a succinct theory of scattering 
and introduce an "inside out" scattering operator of a rever- 
berating environment. 

In section [HI] we show how a change in the environment 
of a linear multi-port system (typically, the interconnect part 
of some electronic system) influences the multi-port model 
parameters of the system. The same analysis can be carried 
through for a multi-port model of Scattering, Thevenin, Nor- 
ton or Hybrid type. Therefore, we shall study a generic inho- 
mogeneous linear model. The essential result of this section is 
that the induced model parameter perturbations are expressed 
in terms of the electromagnetic field scattering operator of that 
environment. It is convenient, but not necessary, to choose 
free space as a reference environment. In that case, we obtain 
expressions for the deviation of the multi-port model param- 
eters due to any object in the system's environment showing 
contrast with respect to free space. The environment scatter- 
ing operator is then defined in the usual way as related to scat- 
tering with respect to vacuum. 

If the environment is considered as a stochastic scatterer, 
the derived relation implies that there are stochastic fluctu- 
ations in the multi-port system's S-parameter model. In or- 
der to characterise these fluctuations, we need a hypothe- 
sis about the nature of the stochastic environment. In sec- 
tion [IV] we give the definition of what we call a Statistically 
Isotropic Environment (SIE). This postulate implies some in- 
teresting properties of the scattering operator of the environ- 
ment. In particular, we find that the variances of the coeffi- 
cients of a matrix representation of the scattering operator sat- 
isfy var(S pil ) = (1 + 8 P q)/N, where N is the finite dimension 
of the wave space in which the fluctuating part of the envi- 
ronment scattering operator works. Observe that this means 
that the scattering operator of the environment satisfies equa- 
tion |[T). At the time of developing our theoretical analysis, we 



were not aware of the fact that such statistically isotropic scat- 
tering problems had already received much attention in the lit- 
erature on neutron scattering. In fact, the formula describing 
the structure of the variance of scattering coefficients is known 
as the Hauser-Feshbach formula . The theory of this formula 
has been formulated in various ways in the past (in particular 
the analysis of Mello and Friedman has many apparent simi- 
larities to our approach 14J]), however, our approach seems to 
lead to the same result in a more direct way. In the context 
of optical scattering in random media, the phenomena of en- 
hanced backscattering leads to similar relations (|<^| Y. Kuga 
and A. Ishimaru, J. Opt. Soc. Am. vol. 1, 1984 p. 831). 
In addition, we use this elementary property of statistically 
isotropic scattering to find variance properties of parameters 
of models which are quite different in nature. 

In section|Vj we show that, under certain circumstances, the 
properties of the environment scattering operator are inherited 
by the fluctuations in the parameters of a generic model for a 
multi-port system in that environment. This section, there- 
fore, provides the theoretical foundation of the relation given 
in equation ([TJ but includes the results, obtained in an entirely 
different way, by the researchers from the University of Mary- 
land. The conclusion of the Maryland analysis that there is an 
essential difference between scattering parameter variance ra- 
tios and impedance variance ratios is not confirmed by our 
analysis. This maybe due to the fact that our hypothesis of 
isotropic scattering environments is incompatible with their 
hypothesis of wave chaotic cavity environment. 

In the appendices, we present some details of the demon- 
strations. In particular, we study the properties of statisti- 
cally isotropic unit vectors in multi-dimensional complex vec- 
tor spaces, which are essential for the theory of scattering by 
isotropic environments. 



H. BASIC RELATIONS 

We shall develop our theoretical analysis using time har- 
monic electromagnetic scattering theory. We want to study the 
modifications in the behaviour of an electronic system due to 
electromagnetic interactions with it's environment. One type 
of such an interaction is related to the electromagnetic fields 
in the environment generated by external sources. Such fields 
induce signals in the system even when the system itself is 
passive. A second type of interaction is related to the electri- 
cal currents, carried by the active system itself, which gener- 
ate electromagnetic fields in the system's environment. These 
fields are scattered by the environment and, when partially re- 
flected back to the system, induce signals which again modify 
the behaviour of the system. The latter situation is, in fact, the 
general one, because we can incorporate the source of some 
ambient field by extending our system to include an additional 
port on which a source is applied. 

We consider configurations, of the general type illustrated 
in Fig. [2] where some electronic system, occupying a domain 
represented by £2 C K 3 , resides in an environment, occupy- 
ing the complementary domain represented by £2 = ffi \ £2 r. 

The "environment domain," £2 + , is the unbounded comple- 
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ment M? \ £2 and contains either a reference distribution of 
reciprocal material or the actual one. In order to simplify the 
presentation of the ideas, we shall take the reference environ- 
ment to be vacuum and suppose that d£l~ has an open neigh- 
bourhood on which the vacuum Maxwell equations hold. 

The "system domain," £2 , is further decomposed in an in- 
terconnect subdomain, in D C £2 , and various port-regions, 
D.P = £2 \D. The complete system consists of, possibly non- 
linear, electronic components, located in £2p, connected to 
the interconnect subsystem. On each connected component 
of £2/>, the low frequency approximation of the electromag- 
netic field is assumed to hold. In the interconnect domain, 
we suppose that all dielectrics and conductors are linear and 
reciprocal. 




FIG. 2: Electronic system decomposed as a linear multi-port inter- 
connect system in D terminated with electronic components in Q.p. 

The electromagnetic fields in the configuration satsify the 
frequency domain Maxwell equations, 

V xH-jaeE = J 

V x E + jcaiiH = 

the constitutive coefficients, e and jj. are supposed to be real 
symmetric tensor functions, representing lossless reciprocal 
media. In the following sections, we concentrate on the model 
of the interconnect system and shall only consider fields corre- 
sponding to port excitations. That means that the distributions 
of electric current source, J, vanish outside £2/>. 

In particular, in the interconnect system domain, D, we use 
these properties to establish the Lorentz field reciprocity rela- 
tion. If {E" ,H a } and {E b ,H b } are two electromagnetic states, 
satisfying the same source-free Maxwell equations in D, we 
have, 

/ v-(E"xH b -E b xH a ) 

JdSlp 

= [ v-(E a xH b -E b xH a ) (2) 
Jdn- 



without any constraints on the nature of the states in the exte- 
rior of D, i.e., the relation holds irrespective of the coefficients 
of the Maxwell equations in £2p and £2 + . In order to simplify 
the notations, we shall write electromagnetics fields as 

y= (E,H) 
and denote forms of the Lorentz type by, 

■&da~ ^\ x V b ) = f v-{E a x H b -E b x H") 

JdQ,- 

We shall use standard scattering theory for the scattering 
description of the environment. This is based on a decom- 
position of the traces on <9£2~ of any electromagnetic field, 
which can exist in the configuration, into two components, 

where V F + is the boundary limit of a solution of the source- 
free vacuum Maxwell equations (and outgoing radiation con- 
dition) in £2 + and is the boundary limit of a solution of the 
source-free vacuum Maxwell equations in £2 . The fact that 
the sum of the two constituents corresponds to a field which 
satisfies a given system of Maxwell equations in £2 + , makes 
that they are related by, 

¥~ = S¥ + 

where S is the scattering operator of the environment. 

Note that, as a consequence of their definitions, we have for 
any pair of states of the same type, 

J? dn - (W a '+ , ) = = J? dn - (^"•- , v*- ) 

In the theoretical development of the following sections, 
we shall encounter Lorentz type integrals, J^sq- , x ¥ b '~), 
where corresponds to a field satisfying the vacuum 

Maxwell equations in Q. and I" 3 to a general system of 
Maxwell equations. This is a frequently occurring case and 
the Lorentz type integral is equivalent to a convenient domain 
integral representation, 

^ n -0P",^)= / (E b -J a ) 
Jn- 

where, assuming a configuration with dielectrics only, J" = 
J+ [jco(e — £q)]E" is called the electric contrast source distri- 
bution of the state 



III. DETERMINISTIC INTERACTION 

In this section, we present the theory of the influence of 
a change in the environment on the model parameters of a 
linear multi-port system. We shall follow the general strategy 
as outlined in ifioll . We represent the linear multi-port system 
by a mathematical model and study the relation between the 
model parameters and the properties of the environment (as 
represented by the model discussed in section HB. 
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The frequency domain model of a linear Alport system has 
the following general form 

y = 77 +Ax 

where x,y, 77 £ and A : C N — ► a linear mapping. The 
vector x, represents the the excitations of the multi-port sys- 
tem and the vector y represents the responses of the multi-port 
system. The inhomogeneity, or "source" term, 77 represents 
the response of the system when there are no excitations on the 
system's ports, i.e., 77 represents sources internal to the multi- 
port system. Note, again, that viewed from the ports of the 
multi-port system the complete environment of the intercon- 
nect system is considered as being part of the multi-port sys- 
tem. This implies that the "internal" source 77 represents elec- 
tromagnetic fields generated by sources in the environment 
coupling to the interconnect system. The Thevenin, Norton, 
Hybrid and Scattering parameter models are all of this form 
and we can develop our theory independently of the specific 
type chosen. If needed, we shall label certain objects with A 
which could then be replaced by Z for a Thevenin model, Y 
for a Norton model H for a Hybrid model or S for a scattering 
parameter model to specialise to a given case. 

The model parameters are 77 and A. These are the basic 
quantities we are interested in. For the purpose of this paper, 
though, we shall only consider the matrix A, because, in the 
context of reverberating environments, all the information we 
need follows from A-parameters. Indeed, if we include the 
input port of some emitting antenna in the environment as an 
additional port of our multi-port system, we can suppose that 
there are no sources in the environment of the resulting ex- 
tended multi-port system, therefore we take 

77 =0 

in the rest of this paper. 

We consider a given A -parameter model of a multi-port sys- 
tem applying to its intended operational environment, say, for 
simplicity, free-space. It is now necessary, to know how a 
deviation from this intended environment translates into de- 
viations of the model parameters. We will first carry out this 
analysis in a deterministic context where the perturbation is 
supposed to be completely known. 

From now on, we shall denote the ideal (reference) model 

by 

y=A°x (3) 

and the true model by 

y = Ax (4) 

A change in an environment can be modeled naturally as an 
electromagnetic scattering problem. So the aim of this section 
is to establish the relation between the scattering model of 
some electromagnetic environment and the matrix A of the 
system's multi-port model. 

We shall now study electromagnetic field states in the con- 
figuration depicted in Fig. [2] Recall that <9£2~ is a closed, pos- 
sibly multi-component, surface, separating the system from 



its environment and d£lp the union of the surfaces enclosing 
the various port regions of the system. The domain D between 
these two surfaces contains the conductors and dielectrics of 
the interconnect system, whereas the domain exterior to the 
surface d£l~ contains all deviations from the intended free 
space environment of the system. 

Upon substitution of the quasi-static approximations on the 
topological components of £2p into the field reciprocity rela- 
tion (O, we obtain (see appendix[A} : 

x bt y a -x at y»=J? da -F> a ,V b ) (5) 

The Lorentz reciprocity relation (0 allows us to find in- 
tegral representations for the difference AA = A — A , if we 
choose the two states appropriately. In our case, 

_ vpi _ vj/0 correS p 0nc j s to x" k = 8^ p a unit excitation ap- 
plied to one of the ports of the system in a free-space 
reference environment, 

- m b = *p a true state, where *P ? corresponds to x\ = 
8k- q , a unit excitation applied to one of the ports of the 
system in its true environment. 

Substituting the model equations, (O and into the equa- 
tion (0, we obtain for the above choices, 

kl 

and using the symmetry of the matrices A (reciprocla media 
outside the port-regions £lp) in the two states 

AA pq = J? ga -(¥ p) V q ) (6) 

Following the scattering theory developed in section H2 the 
field m q can be decomposed on dQr into a part = m®, 
satisfying the same equations as X P® in the exterior of <9£1 , 
and a part satisfying the free-space equations in £l~. In 
addition, there exists a scattering operator between these fields 

Substituting this into © and using Jfg a - ( v F°, l I , °) = 0, we 
get 

AA pq = J? aa -(¥ p ,SW° q ) 
This relation can be discretised using a basis { V^}, 
AA pq = £jS? aa - Wk)Su^ d a- K), Vf ) 

kl 

= ^L pk SkiL q i (7) 

kl 

where we introduced the port-related linear forms L p defined 
by L p (\j/) = JZ'dQ.-^piV)' which have coefficients L pk = 
In-( e k 'Jp) m tne chosen basis {%} and J® is the contrast 
current distribution corresponding to an excitation x k = 8k~ p 
in the free-space environment (Note that in these expressions 
(j,e) = f (j ■ e) has no complex conjugation and should not 
be confused with the L 2 inner product (a,b) = (a,b)). This 
last equation constitutes the searched for relation between the 
perturbation of the A-parameters and the scattering model of 
the environment. 
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IV. SCATTERING OPERATOR OF A STATISTICALLY 
ISOTROPIC REVERBERATING ENVIRONMENT 

In the previous section, we established the relation between 
the perturbations of a linear system's multi-port A-parameter 
model and the scattering model of that system's environment. 
In practice, the true environment of a system is not completely 
known. It is possible to account for this by replacing the scat- 
tering operator of the environment by a stochastic operator. In 
principle, one would like to derive the characteristics of this 
stochastic operator through mathematical analysis of scatter- 
ing problems. However, with the present state of mathematics 
it seems impossible to obtain such properties from first princi- 
ples. Therefore, we shall postulate a minimal set of properties 
of this scattering operator, such as to reflect in the most natural 
way what we understand by a statistically isotropic reverber- 
ating environment. 

In the first place, we identify a decomposition of the en- 
vironment's scattering operator into a fixed average part, S av , 
and a fluctuating part with mean zero, 

s = s av +s flc ' 

The fixed part induces, according to the preceding sections, 
a fixed average deviation from the multi-port model as com- 
pared to the wanted ideal environment. The fluctuating part is 
the one we are interested in here. In analogy with what can be 
found in canonical geometries (see H), the fluctuating part 
of the environment scattering operator is taken to be finite di- 
mensional. This finite dimension grows with frequency in a 
way which depends on the actual geometry of the environ- 
ment. This can be understood physically as being due to the 
fact that in reverberating environments there is a finite number 
of propagating modes reaching the fluctuating parts of the en- 
vironment and coming back into the domain Q. . The remain- 
ing part of the (countably infinite) number of basis functions 
on <9i2~ , does not radiate into the fluctuating environment and 
hence does not play a role in the fluctuating part of the scatter- 
ing operator. An example of this can be found in conventional 
mode stirred chambers, where the mode stirrer is in one part 
of the chamber and only linear combinations of a very specific 
and finite set of plane wave fields, i.e., those corresponding 
to propagating wave guidemodes reflected by a short-circuit 
wall, give non-negligeable fields incident on the mode stirrer. 

In the following development, we shall only talk of this 
zero-mean, fluctuating, part of the total scattering operator of 
the stochastic environment which, for simplicity, we continue 
to write as S. The dimension of the complex propagating wave 
subspace, spanned by the "propagating modes," will be de- 
noted by N. 

Definition IV. 1. A Statistically Isotropic Environment (SIE) 
is such that the fluctuating part of the scattering matrix, rep- 
resenting the contrast of the environment with respect to a free 
space environment, constitutes a spectrally isotropic stochas- 
tic unitary matrix. 

This definition uses a specific stochastic matrix, tailored to our 
needs. 



Definition IV.2 (Spectrally isotropic stochastic matrix). Let 

S : C N — ► C N be a stochastic linear transformation. We call 
this transformation spectrally isotropic if it has a spectral rep- 
resentation 

X 

with stochastic multipliers s% and stochastic rank 1 operators 
of the form Tl^ = v^v^ (transpose with respect to the her- 
mitean inner product), where each v^ is an isotropic unit vec- 
tor ( see appendix\B\). 

Moreover, all the spectral multipliers, s^, and the vectors, 
vx, are mutually statistically independent. In addition, we as- 
sume that all eigenvalues are identically distirbuted stochastic 
variables, i.e., VA E(a^) = and E(\a\ | 2 ) = p 2 - 

In the context of environment scattering, p will be called 
the environment's effective reflection coefficient. 

The inspiration for the above definitions comes from the 
interpretation we can give to a spectral decomposition. A 
linear combination of waves outgoing from the domain £2~, 
which lies in a single eigenspace of the scattering operator, 
will be reflected by the environment as exactly the same lin- 
ear combination modulo only a scalar multiplier s^ (9). If we 
choose the environment to be perfectly isotropic, we would 
require that there exist no preference whatsoever of the vari- 
ous eigenspaces for any subspace of and, moreover, that 
the eigenvalues are statistically the same for any eigenspace. 

An important consequence of this postulate is derived from 
the properties of the cartesian components of isotropic unit 
vectors. 

Proposition IV.3. Let z be a stochastic unit vector, in N di- 
mensional complex vector space, such that the real and imag- 
inary parts constitute a uniformly distributed vector on the 
unit sphere in 2N real dimensions. Then the modules of the 
complex cartesian components, Zk, satisfy: 

E(W 2 ) = i 

E (N 4 )=2^ + ^(AT 3 ) jftf-oo 
E(| Zit | 2 | z ,| 2 ) = E(| Z *| 2 )E(| Z/ | 2 ) + 0(tT 3 ) ifk ^landN - « 

Proof. Let Zk = x + iy, with x and y cartesian components of 
an isotropic unit vector in M? N . We compute, 

E(\z k \ 2 )=E(x 2 +y 2 )=E(x 2 )+E(y 2 ) 
using the results of appendix ?? 

2 _ 1 

~ 2N ~ N 
E(\z k \ 4 )=E(x 4 +y 4 + 2x 2 y 2 ) 

= A r+2E(x 2 y 2 ) 

2JV(2JV + 2) y y ' 
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for the second term on the right hand side, we use the asymp- 
totic estimate from appendixlBlpropositionD 



2 — — + &(N 

Tin*? v 



-3\ 



and, we also have 



2N(2N + 2) 2N2N 
m * m = ^ 7 + {?(N-i), so 

= 2/N 2 + (ff(N- 3 ) 



The third relation can also be shown by straightforward com- 
putation. Let Zk = x + iy and zi — u + iv, then, using proposi- 
tion lB.8l again. we obtain 

E(\z k \ 2 \z l \ 2 )=E((x 2 +y 2 )(u 2 +v 2 )) 

= E(x 2 u 2 +x 2 v 2 +y 2 u 2 +y 2 v 2 ) 

= E(x 2 )E(« 2 ) + ... + <?(-!) 



E(\z k \ 2 )E(\z,\ 2 ) + 



N 3 ' 



□ 



Corollary IV.4. The modules of the complex cartesian com- 
ponents of a complex isotropic vector z € C N satisfy asymp- 
totically 

E(|z,| 4 ) = 2Efc| 2 ) 2 + 0(N~ 3 ) ifN - - 

The scattering matrix model of a SIE satisfies the following 
important relations. 

Proposition IV.5. Let S be the (fluctuating part of the) scat- 
tering matrix of a SIE. Then, 

i) E(S pq ) = 

ii) E(S kl SZ) = (8 k _ m 8,-n + 8k-n8,- m )^ + &{N~ 2 ) 
Proof. As to the average: 

Vp,q E(S w )=£E(v£A*) 
X 

because the average of a sum is the sum of the averages, using 
the statistical independence between s^ and the eigenvectors 
we obtain 

Mp,q E(S m )=£E(v£v*)E(j*) = 
X 

The second property needs a bit more work, 



E{S kl S mn ) = £E(vJ/vf vSjMvff) 
Xfi 

using again the statistical independence between eigenvalues 
and eigenvectors, we find 



E(S kl S mn ) = £E(v£vfv£v£)S^ E(|/ 1 2 ) 



As the joint probability of the cartesian components is a tensor 
product of even functions (see Appendix ??), there are only 
non vanishing contributions on the right hand side when the 
average is over even functions. Therefore, we get 



E(S*/S™)=E 
x 



E(|v* |V | 2 ) S^W 2 

n\^n\M n \ 2 )^-n8l- m p 2 k^m 



using the asymptotics of proposition IIV.3I and 
\\m N ^„Y,x ^(AT 3 ) = G(N~ 2 ), we obtain 



E{S kl S mn ) = (8 k - m 8i- n + 8 k - n 8i- m )^ + 0{N- Z ) 
This yields the proposed relation. 



□ 



In fact, the proof given above shows that the structure of the 
matrix of variances of the scattering coefficients is particularly 
simple in the case of a SIE. 

Corollary IV.6. The variances of the scattering matrix of a 
SIE satisfy 



i) with p ^ q, var(5 P9 ) = \vw(S pp ) 
N — > oo, 



var(S qq ), when 



ii) var(5 M ) - (Sp^ + l)^ + ff{N~ 2 ) 

where N is the dimension of the wave space in which the scat- 
tering matrix is expressed. 

The essential result of this section is that the fluctuating part 
of the scattering operator of a Perfectly Isotropic Reverberat- 
ing Environment (SIE) satisfies the relation (Q3. It is the pur- 
pose of the next section to show that, in a certain way, the 
induced fluctuations of multi-port S-parameters in such an en- 
vironment inherit this basic relation. 



V. INDUCED PERTURBATIONS OF MULTI-PORT 
MODEL PARAMETERS 

With the definition of a SIE given in the preceding section 
and the results of section [TTT1 we are sufficiently equipped to 
characterise the nature of the stochastic perturbations of the 
S-parameter multi-port model of an equipment. In section Hill 
the variations of the model parameters of a linear system, due 
to a variation of the electromagnetic environment, have been 
expressed in terms of bi-linear combinations of scalar prod- 
ucts: 



A/1 



Y^L pk S kl L ql = ^(4:e- k )S kl (J° q: ej) (8) 

kl k .l 



In this expression, {/?} is a set of current distributions, com- 
pletely determined by the system in free space, and {e^ } is a 
set of incident waves, chosen as a base for expansion of regu- 
lar fields in the system domain £2 . 

The variances of the coefficients of the scattering operator 
in a SIE satisfy a simple relation (seeliVli. Because the fluc- 
tuations of the multi-port system's 5-parameters reflect prop- 
erties of the SIE through the latter's scattering operator, we 
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might expect to find similar relations between the variances 
of these perturbations. That this follows indeed depends on 
the following result. 

Proposition V.l. Let L be a matrix of which the rows are or- 
thogonal, i.e., Y.kL p kLqk = 5( p _ 9 )||Lp|| 2 , let S be a spectrally 
isotropic matrix then the coefficients of the matrix A = LSLl 
have variances satisfying the relation 



lim var(A p? ) = ^lhn ^ var(A w _,) var(A w )/4 

where N is the dimension of the subspace in which the fluctu- 
ating part of the scattering operator works. 

Proof. We first show that the averages vanish. Substitution of 
the definitions gives 



l(A pq )=E(£L pk S kl L ql ) 

lc,l 



k,l 




In fact, the orthogonality of the linear forms L p and L q is the 
orthogonality of the system's radiation patterns correspond- 
ing to the excitations of the respective ports p and q. This is 
equivalent to simple power-additivity of the resulting current 
distributions (see also |§]). In other words, when the radia- 
tion power of the simultaneous excitation of ports p and q is 
the sum of the powers of the individual excitations, the pro- 
posed relation should hold (Here, we are neglecting the fact 
that for power additivity only the real part of the inner product 
of radiation patterns needs to vanish!). 

We leave it as an open problem to find practical character- 
isations of situations where this happens. According to the 
variety of configurations where the result has been observed 
experimentally, the conditions are expected to be quite weak. 
This may be due to the fact that the relation is essentially an 
asymptotic relation and that in many configurations N grows 
as / 3 , i.e. the frequency need not be very high before the 
asymptotic estimates hold. 



VI. ESTIMATING THE COUPLING BETWEEN PORTS IN 
ISOTROPIC STOCHASTIC ENVIRONMENTS 



Therefore, the variances are computed as 

var(A w ) = E( L pk SkiL q tL pm S mn L„i) 

k,l,m,n 

= 52 LpkSklLqiLpmLniEfiklSmn) 
k.l.m.n 

using the results of proposition lIV.51 we get 

= H L p i c L q iLp, n L tjn (8k- m 8i-„ + Sjfc_„5/_ m ) — 



k.l.) 



N 



^J[LpkL q iLpkL q i + L p kL q iLpiL q k] — 
k,i yv 

2 2 p2 

[\\ l p\\ \\ L q\\ +Y* L P kL qi L pi L qi<}-jrf 

k.l v 



The second term on the right hand side vanishes due to the 
orthogonality of the rows except when p = q, then we get the 
same result as the first term on the right hand side. This shows, 

var(A )-/lM 2 PJ 2 P 2 V + <1 
Var(A ^-\2||LjV P = q 

This implies, that with p ^ q, 

var(A /)c/ ) 2 = var(A w ,)var(A w )/4 



which is the result we had to prove. 



□ 



The orthogonality between the forms L p , on which depends 
the theoretical proof of equation (OJ, can be obtained for dif- 
ferent reasons. The relation of proposition IV. 1 1 between the 
variances can therefore appear in various different situations. 



In the variance ratios derived in the previous section, the co- 
efficients which determine the actual fluctuations in the model 
coefficients did not appear. In practical situations, however, it 
is frequently necessary to be able estimate the variance of a 
coefficient describing the coupling between two ports. For 
that we need to be able to evaluate the quantities p and ||L„|| 
appearing in the expressions. 

The essential step is to observe that the norms of the port- 
bound linear forms L p give the free-space radiation power of 
the system under unit excitation at its p-th port. We elaborate 
this for the principal multi-port models to which our theory 
applies. 



\L p \\ — ^rad;p 

|L„|| 2 = G, 



rad;p 

[i-\Sp P \ 2 )c 



(Thevenin model) 
(Norton model) 
(S-parameter model) 



where ^ ra d;p is the port's radiation resistance, G ra d; P the port's 
radiation conductance and (1 — \S PP \ 2 )C the port's mismatch 
times a radiation efficiency coefficient C < 1 . 

An isotropic stochastic environment is characterised by a 
single global coefficient p, which we called the environment's 
effective reflection coefficient. We can determine this coeffi- 
cient experimentally, using the S-parameter model for some 
ideal reference two-port system with perfectly adapted loss- 
less antennas, i.e., S pp = and C = 1, we get 

p 2 = var(5i2) 

because in that configuration we have ||Li || 2 = 1 = j|L2|| 2 . It is 
of course possible to do a more realistic computation when the 
mismatch factors and radiation efficiencies of the antennas in 
the reference siutation are known. Of course, the Thevenin or 
Norton model could also be used if the radiation resistance or 
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radiation conductance of the antennas were known. In prac- 
tice, though, the scattering parameters are the most easily ac- 
cessed by measurement. 

To determine the coefficients \\L p \\ 2 for a port on an arbi- 
trary system, port, we can put the system in a conventional 
mode stirred chamber and measure the variance statistic of 
Spp (or R pp or G pp ) over one turn of the mode stirrer. To- 
gether, with the p obtained with the reference setup in the 
actual stochastic environment, we then have the variance of 
the coupling coefficients between any ports in the stochastic 
environment. 



Conclusion 

In this paper, we have developed a theoretical model for the 
fluctuations in the model parameters of linear multi-port sys- 
tems induced by stochastic reverberating environments. This 
theory corroborates an experimentally observed correlation 
between the variances of the reflection coefficients and the 
variances of the transmission coefficients. In a certain way, 



the fluctuations in the reflection coefficients are shown to mea- 
sure the strength of the coupling between a system port and the 
system's environment. This gives a practical, quantitative, ap- 
plication of a familiar idea from antenna reciprocity, putting 
radiation and reception properties into correspondence. On 
the one hand these results can be used to judge on the qual- 
ity of the statistical isotropy of the environment. On the other 
hand, the derived relations can also be applied, in such statis- 
tically isotropic environments, for estimating induced signal 
levels from only ambient field levels and reflection measure- 
ments. 

The theory developed here should be compared to the de- 
velopment presented in ifTHl . The results presented there con- 
flict with ours. Whereas in our analysis all multi-port models 
are handled in the same way, that reference singles out the 
impedance model as a special case. The apparent incompati- 
bility of the results may be due to the fact that our hypothesis 
stating that the stochastic environment is characterised by a 
statistically isotropic scattering matrix appears to be incom- 
patible with the hypothesis of a cavity with wave chaos as 
used in reference IU1I1 . 
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APPENDIX A: LOW FREQUENCY FIELD 
DECOMPOSITIONS RELATED TO ELECTRONIC PORTS 

In this appendix, we present a field decomposition on d£lp 
and derive the consequences for surface integrals of Lorentz 
type. 

As on each simply connected component of D.p we can 
use the low frequency approximation for the electromagnetic 
field, we have 

Vx£ = 0^£ = -V9 
VxH = / 

where we used the Poincare lemma. (If the component should 
not be simply connected, we can eliminate internal exclusions 
to arrive at a global potential.) 

We want to relate Lorentz type integrals to multi-port S- 
parameters. Therefore, we have to proceed intwo steps: first 
relate fields to the usual multi-port quantities, voltage and cur- 
rent, and then field decompositions to the multi-port wave de- 
compositions used in 5-parameter models. 

The first relation is found by subsituting the low frequency 
approximation into the surface integral 

Jn-{E"y.H b ) = -V a -I b 

where V a is the vector of port voltages and I b a vector of port 
currents. This derivation is classical (see ffl [loll ). 

In the vector space, C 2 " of all pairs (V,/), we define the 
projection operators, 
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by 



IF 



i ±m 

±R~ 1 l I 



(Al) 



where R is some positive real constant. 

Let 0^ = V ± /\/2R where is the voltage component in 
the image of FT . Then using 

V = V~+V + 
1 = 1 +I + 



and 7* = ±R l V ± we easily compute 



.y a .i b + y b . i" 



The results of this appendix can then be summarised in the 
following equation 



n ■ (E" xH b — E h x H") 



. 0*;- _ <^;+ . 0« ; - (A2) 



Observe that this relation is the low-frequency analogue of the 
field decomposition used in electromagnetic scattering theory. 



APPENDIX B: SOME RESULTS FROM PROBABILITY 
THEORY 

In this paper, we use classical probability theory as pre- 
sented in textbooks. However, the specific results we need 
our theory are not easily found in the literature. Therefore, we 
feel obliged to present a succinct presentation dedicated to the 
analysis of stochastic unit vectors. 

Mathematical probability theory is formulated in terms 
of Borel measures on topological spaces. In applications, 
though, a probability space is frequently given as a differen- 
tiable n-dimensional manifold, X, and the probability measure 
as a "volume element," i.e., an exterior differential n-form, 
T £ A", on this manifold. If this is the case, we shall speak of 
a probability manifold, written as (X,x). the measure of a set 
U C X is then represented by (t, U) = jjj T. 



Marginal probabilities for charts on probability manifolds 

Charts on a manifold define coordinate functions on open 
subsets of the manifold. Therefore, the individual coordinates 
can be interpreted as stochastic variables (after "normalisa- 
tion" relative to the measure of the open subset in question). 
The probability that a certain coordinate will be in a given 
interval of K is the probability measure of the subset of the 
manifold for which the points have the chosen coordinate in 
that interval. 

Definition B.l. Let <p : X D U — > V c W be a chart on a 
probability manifold (X,t) of dimension n. The probability 
that (p k e A C R is defined by 

F q>k (A) = (T,( P - 1 [Vn(n k )- 1 [A}})/(T,U) 



where 7Ij :W 3 xt-^r x k is the projection on the k-th component 
(note, (p k = K k o (p and hence we might have written ((p k )~ l , 
but we prefer the explicit representation here to make appear 
the intermediate subset V H (n k )~ l [A]). 

We can establish a marginal probability density in the fol- 
lowing way 

Proposition B.2. The k-th marginal probability density, jj,^ € 
A", of a chart, <p : X D U — > V C W, on an n-dimensional 
probability manifold, {X,x), is defined by 

VA c R (m£,A> b = (t, (p~' [V n 7Zf7 ' [A]])x/ (t, U)x 

Proof. In fact, the only new thing with respect to the above 
definition is the fact that densities are uniquely defined by 
their evaluation on integration domains. □ 



Isotropic stochastic unit vectors 

In this section, we investigate the properties of what we 
call isotropic stochastic unit vectors in R". (This is not a real 
restriction as a unit vector z in Hermitean C", with ||z||c n = 
VEk^izk) 2 +3(z/t) 2 , is also a unit vector in Euclidean R 2 ".) 

Definition B.3. A stochastic unit vector in W is defined by 
a probability manifold (S"~~ ,t), where (TjS" -1 ) = 1 and T 
some density on the unit sphere S"~ l . An isotropic stochastic 
unit vector is defined by the normalised measure on the unit 
sphere induced by the, rotation invariant, Euclidean norm on 
W, t c „-i, i.e., t = t*-i/|S7 _1 |. 



Proposition B.4. The cartesian components of an isotropic 
stochastic unit vector in R" are stochastic variables z € 
[— 1,+1] with a probability density p{z)dz, and 

P (z) = (i -z 2 )("- 3) / 2 r(»/2)/V^r((«- 1)/2) 



Proof. We use the expression of proposition IB. 21 for the 
marginal probability measure of probability measures given 
as densities on manifolds. We obtain, VA C (—1,1), 



»— 2j 

^rndz,Ah 



-n-3 |5"- 21 



■z 



\S" 



Lg = r(n/2)/Vir((n - l)/2), we get, 



((1 -z 2 )("- 3 V 2 r(«/2)/^r((« - l)/2)dzAh 



□ 



We obtain the first few moments 
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Corollary B.5. A cartesian component, x k , of an isotropic 
unit vector in x £ R" has the following moments, 

m\(x k ) =0 
. , 1 

m 2\ x k) = ~~ (equi-partition of variances) 
n 

m 3 (x k )=Q 

i 

m 4 (x k ) = 



n(n + 2) 



Proof. The moments have the following integral representa- 
tion, 



h=-\ sj%T{{n-\)l?) 



As the integrand is an even function on (—1,1), the odd mo- 
ments all vanish. The even moments can be found as twice 
the integral over (0, 1), which leads to standard integrals. 



= 2 



T(n/2) 



z=o ' 0fT((n-l)/2) 
0rT((n-l)/2) r(n/2) 



2r((« + 2)/2) 0fr((n-l)/2) 
r(n/2) _ r(n/2) 



2T(n/2 + l) «T(n/2) 
1 

The fourth moment is computed in a similar way using 
■i • 30rr((n-l)/2) 



2 I z\l-z 2 )^' 2 dz = 

Jz=0 



4r((n+4)/2) 



using r((n + 4)/2) = |(1 + l)r(f) 



3y^T((w-l)/2) 
n(n + 2)r(n/2) 



□ 



Proposition B.6. The cartesian components of an isotropic 
unit vector in R" tend to normal gaussian variables, 
N(0, -33), when n — ► 00. 

Proof. With the following two standard results 



2^-3)/2 -; 2 (^) 



Hm (l-z z ) 

«— 3 



r((«-3)/2 + 3/2) = fn-3 
r((n-3)/2 + l) V 2 

2 



we immediately get 



n^3 
2 



lim (l- z 2 )(n -3)/2 r^/2 



v^r((«-i)/2) V2H 



<7\/27r 



e 2CT 2 



witha 2 = l/(«-3). □ 
Proposition B.7. 77ie joint probability distribution of any m 
cartesian components of an isotropic stochastic unit vector in 
n real dimensions, is given by the probability density function 

I an— m— 1 1 

P (x l , . . . ,x m ) = ^-^n-- 1 (1 - (x «-*)2)<»-*- 3 >/ 2 

Proof. The proof follows the same reasoning as for the ID 
marginal distribution. We start with the joint probability of 
two cartesian components. We know the density on the unit 
sphere S n ,~ , in natural coordinates, 



T^-i = Tw,- 2 Ad9" 1 

A l \in(9"-') 



substitution of the corresponding expression for the first factor 
gives 



sin(0"- 1 )"- 2 T,„- 3 Ad9"- 2 Ade n ~ l 



sm(e"- 2 ) 
in— l\n— 2 „:„/nn—2\n—3. 



-2 A joB-1 



sm(G n - l ) n - z sm(9 n - z ) n - i T s „-3 Ad9"- Z Ad9 



Substitution of a" = cos (9"- L ) then yields 

= (1 - (x") 2 ) ( "" 3)/2 (l - (x"- 1 ) 2 )*"- 4 '/ 2 

Adx"- 1 Adx" (Bl) 

Using this result in the expression for marginal probability 
densities, we get, for any two cartesian components (The lack 
of symmetry in the joint probability density function is related 
to the fact that the components of an isotropic unit vector are 
not statistically independent and their squares have to sum up 
to unity.) 

/7(Jc n -^x")^/x n - 1 ^ ix" = (l-(x") 2 )("- 3 )/ 2 (l-(x"- 1 ) 2 )(' , - 4 )/ 2 



rfti— 3 1 



r.n-1 



dx"- l dx n (B2) 



It should be clear from the above procedure, that we can do 
this for any number of dimensions and get the expression of 
the proposition. □ 

Proposition B.8. Let x and y be two different cartesian com- 
ponents of an isotropic stochastic unit vector in R", 

E(xy) = 
E(x 2 y 2 ) = E{x 2 )E(y 2 ) + e{n^) 

Proof. The first relation is a simple consequence of the sym- 
metry properties of the joint probability distribution. The 
complete proof of the asymptotic estimate in the second re- 
lation is very technical and too long for this paper. Here we 
content ourselves by observing that the limit distributions are 
Gaussian and, hence, the proposition gives the correct limits. 

□ 



